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Abstract 

The bond dissociation energies ( D e ) of C2H2, C2H4, C2He, N2, N2H2, and 
N2H4 are studied at various levels of correlation treatment. The convergence of 
D e with respect to the one-particle basis is studied at the single-reference modified 
coupled-pair functional (MCPF) level. At all levels of correlation treatment, the er- 
rors in the bond dissociation energies increase with the degree of multiple bond char- 
acter. The multireference configuration-interaction (MRCI) D e values, corrected for 
an estimate of higher excitations, are in excellent agreement with those determined 
using the size-extensive averaged coupled-pair functional (ACPF) method. We find 
that the full-valence complete-active-space self-consistent-field (CASSCF)/MRCI 
calculations are reproduced very well by MRCI calculations based on a CASSCF 
calculation that includes in the active space only those electrons involved in the 
C-C or N-N bonds. To achieve chemical accuracy (1 kcal/mole) for the D e val- 
ues of the doubly bonded species C2H4 and N2H2 requires one particle basis sets 
including up through h angular momentum functions (Z= 5 ) and a multireference 
treatment of electron correlation: still higher levels of calculation are required to 
achieve chemical accuracy for the triply bonded species C2H2 and N2. 
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I. Introduction 

The modelling of combustion processes requires accurate bond energies 
for many molecular species. Recent work [1-3] has shown that C— C and C— H 
single bond dissociation energies can be straightforwardly computed using current 
ab initio methods. However, for electron-dense systems such as C 2 H 2 and N 2 , 
it is well known that extensive one-particle basis sets and correlation treatments 
are required to approach chemical accuracy for the binding energy. The greater 
degree of difficulty in treating multiply-bonded systems can potentially complicate 
the determination of heats of formation. For example the Haber process for the 
synthesis of ammonia: 

N 2 + 3H 2 2NH 3 (1) 

involves breaking a triple bond and three H-H single bonds and forming six polar 
N-H single bonds. Thus if the triple bond is treated less accurately than the single 
bonds, this can potentially introduce an error into the calculated heat of formation. 
Yet calculations using quadratic configuration-interaction methods in conjunction 
with extrapolation to the complete basis set limit (CBS-QCI) [4] give excellent 
agreement with experiment [5]. On the other hand, the computation of even C- 
H single bond dissociation energies is complicated if in the process of breaking 
the bond the degree of multiple bond character elsewhere in the molecule changes 
significantly [3]. Thus it is important to understand the limitations of various 
theoretical methods as to their ability to treat equivalently single, double, and 
triple bonds. 

In this work we determine the binding energies of the singly bonded species 
C 2 H 6 and N 2 H 4 , the doubly bonded species C 2 H 4 and N 2 H 2 , and the triply bonded 
species C 2 H 2 and N 2 at various levels of theory. One goal of our study is to calibrate 
the single-reference-based MCPF approach, as this can be applied to relatively 
large systems using high quality one-particle basis sets. For our multireference 
calculations, we consider whether it is necessary or desirable to perform full valence 
CASSCF calculations, or whether use of a smaller active space including only those 
electrons affected in the bond-breaking process is adequate. In the next section 
we discuss the details of the theoretical approaches employed in this work. In 
section III we compare the calculated binding energies at various levels of theory 
with experiment, as accurate Z?o values and zero-point corrections are available for 
these species. Our conclusions are presented in section IV. 
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II. Methods 

Except for basis set calibration studies, we have used atomic natural or- 
bital (ANO) basis sets [6] of the form (135 8p 6<f 4/)/ [45 3 p 2 d If] for carbon and 
nitrogen and (85 6p 4d)/[4s 2 p Id] for hydrogen. These basis sets axe described in 
more detail in Refs. 6 and 7. For all systems we carried out MCPF calculations 
after adding an additional 5, p, d and / ANO function, as well as a { 2 g)/\ig] ANO 
function to the heavy atom basis. In addition, for the doubly and triply bonded 
species we added an additional 5 — g ANO and a primitive h function. The orbital 
exponents of the h function added to C and N are 1.24 and 1.728, respectively. Thus 
the largest contracted basis sets employed for C and N were [65 5 p 4 d 3/ 2 g lh]. 
We also expanded the hydrogen basis to a nearly complete [65 5 p 2 d If] set [8] 
for calculations on C2H4. Extehsive basis set studies were also considered for N2 
in earlier work [9,10]. The smaller [45 3p 2 d If] basis sets, however, are entirely 
adequate for addressing the limitations of the various correlation approaches. Only 
the pure spherical harmonic components of the basis functions were used in the 
calculations. 

Both single- and multireference-based correlation approaches were used to 
determine the bond dissociation energies. For our single-reference-based approach, 
we use the approximately size-extensive MCPF method [11], which has been shown 
to provide accurate dissociation energies for single bonds provided both the molecule 
and its fragments are well described by an SCF reference. This is the case for all of 
the systems considered in this work. 

Our multireference-based methods include the complete-active-space self- 
consistent-field (CASSCF) multireference configuration-interaction (MRCI) and av- 
eraged coupled-pair functional (ACPF) methods [12]. In the initial choice of the 
active space we included all of the bonding and antibonding orbitals and electrons. 
Thus all electrons were correlated except the nominally Is electrons on C and Is 
and 2s electrons on N. However, for C2H6 and N2H4 this resulted in a very large 
CASSCF-CI expansion. For example, for C2H6 the full-valence calculation corre- 
sponds to 14 electrons in 14 orbitals. As this was very time consuming we carried 
out a 10 in 10 calculation instead, where one orbital in each of the a\ g and a 2 U 
symmetries of the Dzd point group were moved to the inactive space and the corre- 
sponding antibonding orbitals were moved to the secondary space. These orbitals 
correspond to the symmetric and antisymmetric (with respect to the C-C bond) 
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combinations of the CH bond orbitals that axe totally symmetric on each fragment. 
We believe that MRCI calculations based on either a 14 in 14 or 10 in 10 CASSCF 
treatment will give very similar results. An analogous reduction of the N2H4 cal- 
culation from 10 in 10 to 6 in 6 was also performed. The second choice of the 
CASSCF active space included only those electrons involved in either the C-C or 
N-N bonds. This was possible as these orbitals could be clearly distinguished from 
those involved in the bonds to hydrogen. For the MRCI and ACPF calculations it 
was not possible to include all single and double excitations from the full-valence 
CASSCF wave function, as this led to prohibitively long expansions. Thus the ref- 
erence lists were restricted to include all occupations for which the absolute value 
of the coefficient of any one of its component spin couplings exceeded a designated 
threshold (generally 0.05) in the CASSCF wave function. An advantage of the 
smaller CASSCF active spaces is that no reference selection was required. That is, 
all MRCI and ACPF calculations involving the smaller active spaces include the full 
CASSCF as the reference wave function. All electrons, except the Is electrons on 
carbon and nitrogen, were correlated in the MCPF, MRCI and ACPF calculations. 
The effect of higher excitations in the MRCI calculations were estimated using the 
multireference analog of the Davidson correction, denoted +Q. The reference wave 
function for the ACPF calculations was expanded until it included all configurations 
with coefficients greater than 0.05 in the final ACPF wave function. 

The D e values were computed as the energy difference between the equilib- 
rium geometry and a dissociated supermolecule structure. Experimental geometries 
were used for both the molecules and fragments [13-15]. For the MCPF treatment of 
C2H2, it was necessary to dissociate to two CH fragments in the excited 4 E” state, 
and then to correct to ground state fragments using the known 4 E“ ^- 2 H separation 
for CH [13]. For N 2 H 2 (diimide), we used the experimental geometry for the trans 
structure. For N2H4 we used the gauche equilibrium structure [14,15]. All calcu- 
lations were carried out using the MOLECULE-SWEDEN [16] program system on 
the NASA Ames Research Center Central Computing Facility Y-MP/832. 

III. Results and Discussion 

The D e values as a function of correlation treatment are summarized in 
Table I for the six molecular systems considered in this work. An important con- 
clusion that can be drawn from the results in Table I is the excellent agreement of the 
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MRCI+Q and ACPF results based on both the full- valence and smaller CASSCF 
active spaces. This is particularly true of the C2H2T1 species where the agreement is 
nearly exact. In fact, there is more variation of the full- valence MRCI results with 
selection threshold than between the MRCI calculations based on the two choices of 
active space. The excellent agreement between the MRCI+Q and ACPF D e values 
indicates that the contribution from higher excitations has been accounted for sat- 
isfactorily. Probably our best results axe the ACPF values using the smaller active 
space, as the ACPF incorporates the effects of higher excitations in a more rigorous 
manner and the smaller active space avoids errors inherent in reference selection. In 
comparing these ACPF values with the MCPF values, we see that the discrepancy 
increases with the degree of multiple bond character. For single bonds the error is 
between 0-1 kcal/mole, for double bonds it is of the order 3-4 kcal/mole, and for 
triple bonds it is of the order of 6 kcal/mole. This is consistent with expectations 
that higher excitations are more important in describing multiply bonded systems. 

Our theoretical results in Table I are of course limited in accuracy by ap- 
proximations in both the one- and n-particle expansions. Previous studies [10,17] 
for N2 indicate that the ACPF method gives binding energies within 1-2 kcal/mole 
of the full configuration-interaction (FCI) limit. We expect that the ACPF treat- 
ment will be even closer to the FCI limit for the singly and doubly bonded species, 
since higher excitations will be less important than for N2. However, the use of a 
[4 s 3 p 2d If] ANO basis set for C and N results in a substantially larger underesti- 
mation of the D c value than does the use of the ACPF method for accounting for 
electron correlation. To assess these limitations, we have carried out basis set cali- 
bration studies at the MCPF level — see Table II. Two trends are apparent from the 
results in Table II. First, the magnitude of the basis set corrections increase with the 
degree of multiple bondedness, and second, the corrections are significantly larger 
for the nitrogen species. It is interesting to note that the basis set requirements 
are also larger for N-H as compared with C-H bonds. For example, at the MCPF 
level, increasing the basis set from [4 s 3 p 2d If] to [5s 4 p Zd 2/ 1 g] increases the 
first N-H bond dissociation energy of NH3 by 1.2 kcal/mole, whereas it increases 
the C-H bond dissociation energy of CH4 by only 0.5 kcal/mole. 

As can be seen from Table II, the effect of expanding the basis from 
[5s 4 p Zd 2/ 1#] to [6s 5 p 4 d 3/ 2 g 1A] is on the average about 0.35 that of expanding 
it from [4s 3 p 2d If] to [5s 4p 3d 2/ 1 g]. Also, the contribution from adding a single 
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h function is about 0.30-0.35 that of adding the first ANO g function. We have made 
an estimate for the remaining basis set incompleteness by assuming that further ba- 
sis set saturation and the addition of still higher angular momentum basis functions 
increases D e in a similar manner. This is essentially a geometric aprroximation: 
the correlation energy itself converges more slowly than this, but the convergence 
of the binding energy is uncertain. Our approximation will certainly be an under- 
estimate of the remaining basis set correction, especially since it does not account 
for further saturation in the primitive basis: note that expanding the N basis from 
(135 8p 6d 4/ 2 g 1 h)/[6s 5 p 4 d 3/ 2 g Ih] to (185 13p 6d 5/ 4 g Zh)/[6s 5 p 4 d 3/ 2 g l/i] 
increases the D e of N2 by 0.2 kcal/mole at the MRCI level [10]. The effect of basis 
set saturation in N2 is very similar at the MCPF and MRCI levels of correlation 
treatment, so that the MCPF approach is expected to be adequate for evaluating 
basis sets effects on the D e of the other molecules considered in this work. Our 
extrapolation procedure also does not account for improvements in the hydrogen 
basis. However, expanding the hydrogen basis set from [45 2 p Id] to [65 5p 2d 1/], 
in conjunction with a [5s 4p 3d 2/ 1 g] ANO basis set for carbon, increased the D e 
of C2H4 by only 0.1 kcal/mole. This very small increase in D e is not surprising, 
of course, since the hydrogen basis functions contribute very little to the C-C bond 
breaking process. The values in Table II in the line denoted “total” are our estimate 
of expanding the heavy atom basis set from [4s 3p 2d 1/] to the basis set limit at the 
MCPF level of correlation treatment. We use these values in Table III to correct 
our ACPF D e values in the [4s 3p 2d 1/] basis before comparing with experiment 
and the values deduced using the G1 approach [18]. We have not included basis 
set superposition error (BSSE) corrections, as these are expected to be small. For 
example, the BSSE determined using the counterpoise correction with the full ghost 
basis was 0.28 kcal/mole for N2 in the [65 5p 4d 3/ 2 g 1 h) basis [10]. This correction 
is expected to be less for the carbon species and for the singly and doubly bonded 
nitrogen species. 

To facilitate comparison with experiment, we have converted our D e values 
to Dq values using the experimental vibrational frequencies [19,20]. The Do values 
denoted theory in Table III are our computed D e values at the ACPF level (using 
the smaller active space) in the [4s 3p 2d l//4s 2p Id] ANO basis corrected for zero- 
point and by our estimate for extrapolation to the basis set limit. Note that we have 
deduced the experimental binding energy of N2H2 using a heat of formation, A2?^ 98 , 
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for trans-diimide of 36±2 kcal/mole [21]. Our calculated Dq value is inconsistent 
with the two larger experimental values [22,23] for AH* 9 * } as well as the calculated 
value of Casewit and Goddard [24]. 

Our theoretical Dq values in Table III tend to be smaller than either exper- 
iment [13,19,21,25] or the values determined using the G1 approach [18]. Overall 
our Dq values for the carbon species are in better agreement with experiment than 
are the Gl values, whereas the opposite is true for the nitrogen species. Our values 
are almost certainly lower bounds, as the ACPF method tends to underestimate 
slightly (at least for the multiply bonded species) the correlation contribution to 
He, and our basis set extrapolation procedure does not account for saturation of 
the primitive basis. 

Considering that our theoretical values in Table III are probable lower 
bounds, the correct Dq value for C2H4 probably lies in the upper half of the experi- 
mental error bars. Most of the remaining 2 kcal/mole discrepancy with experiment 
for the C-C bond dissociation energy of C2H2 can probably be attributed to lim- 
itations in the ACPF method. Our theoretical Dq value of N2H4 is smaller than 
the experimental value, but within the error bars that are relatively large due to 
the uncertainty in the heat of formation of NH2. The Gl value is also smaller 
than experiment, but it should be noted that the Gl value is determined using 
a theoretically determined zero-point correction of —9.5 kcal/mole, instead of the 
experimental value of —7.2 kcal/mole. Thus the Gl Ho value would be larger than 
experiment, in analogy with C2H6, if the experimental zero-point correction [20] 
was used instead. 

Petersson and Mantzaris [4] have computed the heat of formation of am- 
monia using the CBS-QCI model. Their value of — 8.83±0.5 kcal/mole is in good 
agreement with the accurate experimental value [5] of — 9.34±0.01 kcal/mole. They 
present a careful analysis of why the value of —6.52 kcal/mole obtained with the 
generally very reliable Gl model [26] is in error by nearly 3 kcal/mole. Since the 
Haber process for the synthesis of ammonia involves breaking a triple bond in N2 
and three single H-H bonds and subsequently forming polar single bonds in NH3, 
we considered how reliably the MCPF approach could determine the heat of forma- 
tion using the [4s 3p 2d l//4s 2p Id] ANO basis. Using Petersson and Mantzans’s 
value of 10.31 kcal/mole for the zero-point correction, we compute a AH 9 value 
of —6.97 and —8.25 kcal/mole at the SCF and MCPF levels. For comparison in 
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this basis set the CCSD and CCSD(T) methods produce AJI^ values of — 8.83 and 
— 7.81 kcal/mole, respectively. Thus the MCPF value lies between the CCSD and 
CCSD(T) values. Triple excitations tend to decrease A Hj, because they are more 
important for N2. Expanding the basis set to [53 4 p 3 d 2/ ly/43 2p Id], decreases 
the MCPF value to — 9.24 kcal/mole, in neax perfect agreement with experiment. 
Thus there is a relatively small effect of basis set saturation, which in large part 
arises from a cancellation of errors that occurs when the process is formulated as an 
isogyric reaction. As there is very little error in the bond energy of H2, obtaining 
the correct heat of reaction for ( 1 ) (and thus the correct heat of formation for NH3) 
relies on a cancellation of the error in the triple bond of N2 and the error in six N-H 
bonds. This cancellation is nearly exact at the MCPF level in the [ 5 s 4 p 3 d 2/ 1#] 
basis, and based on our basis set studies we expect that the A Hf value of NH3 will 
also be accurate at the MCPF level in the basis set limit. This is consistent with 
the rather good result obtained by Petersson and Mantzans [ 4 ] with the CBS-QCI 
approach. We should also note that it should be possible to determine Ajffy for 
CH4 accurately by a similar procedure. The D e of C2 in the [ 4 s 3 p 2d If] basis at 
the MCPF level is about 11 kcal/mole less than experiment [ 27 ], and the C-H bond 
strength at this level is 1.5 kcal/mole less than experiment [ 25 ]. Thus the errors are 
comparable on both sides of the reaction 4H2 + C2 2CH4. 

IV. Conclusions 

We have carried out calculations for the singly bonded species, C2H6 and 
N2H4, the doubly bonded species, C2H4 and N2H2, and the triply bonded species 
C2H2 and N2 at various levels of correlation treatment to assess how the errors in 
the binding energies vary with the degree of multiple bondedness. We find that 
the errors in the single-reference approach as compared with the multireference 
ACPF approach increase from 0-1 kcal/mole for single bonds to 3-4 kcal/mole for 
double bonds to approximately 6 kcal/mole for triple bonds. In addition, the one- 
particle basis set requirements also increase substantially with the degree of multiple 
bondedness. Thus while it may be possible to compute binding energies to chemical 
accuracy by direct ab initio methods for singly bonded species, it will be necessary 
to either employ very large one-particle basis sets (including at least h functions) or 
an extrapolation procedure to account for the remaining basis set incompleteness, 
such as those performed in the G1 (and G2) approaches or the CBS-QCI model, 
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in order to obtain quantitative results for the D e values of multiple bonds. The 
multireference calculations performed in this work underscore the desirability of 
restricting the CASSCF space to include only those orbitals and electrons involved 
in the bond breaking process. 
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Table I. D e values (kcal/mole) as a function of the level of correlation treatment. 


System 

Active space 



D e 





CASSCF 

MRCI 

MRCI+Q 

ACPF 

MCPF b 

c 2 h 2 

10x10 

232.0 

229.1 

228.5 

228.3 



6x6 

236.4 

227.8 

228.8 

228.3 

223.0 

c 2 h 4 

12x12 

167.2 

175.8 

177.6 

177.4 



4x4 


176.3 

177.5 

177.3 

174.6 

c 2 h 6 

10x10° 


92.2 

94.4 

94.6 



2x2 

80.8 

93.0 

94.4 

94.4 

94.2 

n 2 

6x6 

203.2 

215.7 

215.8 

215.7 

209.3 

n 2 h 2 

8x8 

114.3 

122.3 

123.9 

124.0 



4x4 


122.0 

122.8 

122.8 

118.3 

n 2 h 4 

6x6° 

64.8 

67.2 

68.3 

68.3 



2x2 

51.8 

66.3 

67.6 

67.5 

66.8 


° See the text. , 

6 Based on SCF orbitals. 
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Table II. MCPF D e values (kcal/mole) as a function of the heavy atom basis set. 


Basis 0 


• 


c 2 h 2 

c 2 h 4 


[4s 3p 2d If] 

223.03 

174.56 

[5s 4p 3d 2 /] 

225.35 

175.62 

[5s 4p 3d 2/ lg] 

227.52 

176.76' 

[6s 5p 4d 3/ 2 g] 

228.26 

177.24 

[6s 5p 4d 3/ 2g lh] 

228.88 

177.60 

oo d 

229.47 

178.12 

Total 

6.44 

3.56 


Molecule 


C 2 H 6 '’ 

n 2 

N 2 H 2 

n 2 h 4 * 

94.15 

209.26 

118.33 

66.82 


212.95 

120.15 


95.16 

215.62 

122.04 

69.05 


217.11 

122.77 



217.91 

123.34 


95.8 

219.20 

124.04 

70.3 

1.6 

9.94 

5.71 

3.5 


° The hydrogen basis is [4s 2p Id]. 

* The total basis set correction is based on the difference between the D e values in 
the [4s 3 p 2d If] and [5s 4 p 3 d 2f lg] basis sets and the total correction for the 
corresponding doubly bonded system. 

c The D e value increases to 176.85 when the hydrogen basis set is expanded to 
[6s5p2dl/]. 

d Estimate based on a geometric extrapolation. 
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Table III. Comparison of the calculated C-C and N-N bond strengths (kcal/mole) 
with experiment. 


System 


Corrections 



D o 




Zero-point Basis set* 

Theory 

Gl c 

Expt. 

C 2 H 2 

228.3 

-8.1 

+6.4 

226.6 

227.0 

228.8+0.7* 

c 2 h 4 

177.3 

-9.7 

+3.6 

171.3 

175.7 

171.0+1.2* 

c 2 h 6 

94.4 

-9.0 

+1.6 

87.0 

88.7 

87.8+0.3* 

n 2 

215.7 

-3.4 

+9.9 

222.2 

224.8 

225.0' 

n 2 h 2 

122.8 

-7.8 

+5.7 

120.7 

. . . 

121.8+2.0® 

n 2 h 4 

67.5 

-7.2 

+3.5 

63.8 

64.6 

65.4+2.0* 


° ACPF results obtained using the [4s 3 p 2d 1/ /4s 2 p Id] ANO basis sets. 
h Estimate for basis set incompleteness — see Table II. 
c Reference 18. 

* Reference 25. 

e Reference 19. The error bars reflect the uncertainty in the heat of formation of 

CH 3 . 

f Reference 13. 

9 Reference 21. 

Reference 19. The error bars reflect the uncertainty in the heat of formation of 
NH 2 . 
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